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ABSTRACT 

Aims. We revisit the derivation of the mass function and the bias of dark matter halos for non-Gaussian initial conditions. 
Methods. We use a steepest-descent approach to point out that exact results can be obtained for the large-mass tail of 
the halo mass function and the two-point correlation of massive halos. Focusing on primordial non-Gaussianity of the 
local type, we check that these results agree with numerical simulations. 

Results. The high-mass cutoff of the halo mass function is of the same form as the one obtained from the Press- 
Schechter formalism, but with a linear threshold 5l that depends on the definition of the halo (i.e. Sl — 1.59 for a 
nonlinear density contrast of 200) . We show that a simple formula, that obeys this large-mass asymptotic and uses the 
fit obtained for Gaussian initial conditions, matches numerical simulations while keeping the mass function normalized 
to unity. Next, deriving the real-space halo two-point correlation, in the spirit of Kaiser (1984), and taking a Fourier 
transform, we obtain a good agreement with simulations for the correction to the halo bias, A6Af (fc, /nl), due to 
primordial non-Gaussianity, without any free parameter. The nonlinear real-space expression can be useful to check 
that the "linearized" bias is a valid approximation. Moreover, it clearly shows how the baryon acoustic oscillation at 
~ 100/i~^Mpc is amplified by the bias of massive halos and modified by primordial non-Gaussianity. At smaller scales, 
30 < a; < 90/i~^Mpc, the correction to the real-space bias roughly scales as /nl 6m(/nl = 0) x^. 

Key words, gravitation; cosmology: theory - large-scale structure of Universe 



1. Introduction 

,, Standard single-field slow-roll inflationary models predict 
' a nearly scale-invariant and Gaussian spectrum of primor- 
I dial curvature fluctuations (e.g., Bartolo et al. 2004). This 
. agrees with current observations of the cosmic microwave 
' background (CMB) anisotropies (Komatsu et al. 2008) and 
, of large-scale structures (Slosar et al. 2008). Nevertheless, 
' several inflationary models predict a potentially observable 
' level of non-Gaussianity (e.g., Bartolo et al. 2004 for a re- 
, view), so that constraining or detecting primordial non- 
• Gaussianity is an important task for current cosmological 
studies. This would allow one to rule out some of the many 
inflationary models that have already been proposed. In 
particular, in many cases the non-Gaussianity is of the lo- 
cal type, meaning that it only depends on the local value 
of Bardeen's potential <!>. That is, the latter can be decom- 
posed as 

local type: $(x) - 0(x) + /nl (^(x)^ - (02)), (i) 

where (jj is a Gaussian random field. Simple slow-roll infla- 
tion gives a parameter /nl of order lO^'^ but this would 
be masked by the nonlinearities that arise from the dy- 
namics (e.g., from the nonlinearity of Einstein's equations, 
see Bartolo et al. 2004) or from the physical processes in- 
volved by the observables (e.g., perturbations at recombi- 
nation that affect the CMB, see Senatore et al. 2008), which 
give an effective /nl of order unity. Large values of /nl can 
be obtained for instance from multifield infiation (Bartolo 
et al. 2002; Lyth et al. 2003), self-interactions (Falk et al. 
1993), tachyonic preheating in hybrid infiation (Barnaby & 



Cline 2006) or ghost infiation (Arkani-Hamed et al. 2004). 
Current limits are — 9 < /nl < 111 from CMB (Komatsu et 
al. 2008) and —29 < /nl < 70 from large-scale structures 
(Slosar et al. 2008). 

The effects of primordial non-Gaussianity on large-scale 
structures can be seen for instance through the mass func- 
tion of virialized halos, especially in the large-mass tail as 
the steep falloff magnifies the sensitivity to initial condi- 
tions (Lucchin & Matarrese 1988; Colafrancesco et al. 1989; 
Grossi et al. 2007; Maggiore & Riotto 2009). This allows to 
use the X-ray luminosity function of clusters to constrain 
the amount of non-Gaussianity (Amara & Refregier 2004). 

A second probe of non-Gaussianity is provided by the 
clustering of these halos, measured through their many- 
body correlations. In particular, the halo two-point corre- 
lation can be significantly increased if the underlying pri- 
mordial density field is non-Gaussian (Grinstein & Wise 
1986). More specifically, Dalai et al. (2008) have recently 
shown that primordial non-Gaussianity of the local type ([1]) 
gives rises to a strongly scale dependent bias at large scales, 
whereas in the Gaussian case the bias is roughly constant 
in this range. Thus, at linear order over /nl they obtain in 
Fourier space a correction of the form 

A5M(fc,/NL) = /nl [bM{k,0) 1] J^'^^l]^ly (2) 

where bM{k,0) is the Gaussian-case bias (defined as the 
ratio of the halo and matter power spectra, &|f(fc) = 
Puik)/ P{k), for objects of mass M). Here is the lin- 
ear matter density contrast associated with virialized ob- 
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jects (usually taken as 5^ — 1.686), T{k) is the transfer 
function and D{z) is the linear growth factor, normalized 
as D{z) ^ (1 + z)~^ at high redshift. Then, Dalai et al. 
(2008) checked in numerical simulations that, in agreement 
with Eq. , the halo bias correction roughly grows as l/k"^ 
at low k. This gives rise to a significant and specific signal 
that has already been used to constrain /nl (Slosar et al. 
2008). 

In this article, following a previous work devoted to 
the Gaussian case (Valageas 2009b), we revisit the deriva- 
tions of the halo mass function and of the bias for primor- 
dial non-Gaussianity. Although we focus on the local type 
((1]), our approach also applies to any non-Gaussian model 
where Bardeen's potential can be written as the sum of lin- 
ear and quadratic terms over an auxiliary Gaussian field, 
that is, where /nl becomes a convolution kernel (it also ex- 
tends to cases that contain higher order terms and multiple 
Gaussian fields). 

After introducing our notations and the quantities 
needed for our calculations in section [21 we consider the 
halo mass function in section [31 Here, our aim is to ar- 
gue that the exponential cutoff of the large-mass tail can 
be obtained exactly from a saddle-point approach. This is 
equivalent to the saddle-point computation of Matarrese et 
al. (2000), that is often used to model the non-Gaussian 
halo mass function. However, using a different treatment, 
we simultaneously derive the linear density profile of this 
saddle-point, which allows us to check that the latter is 
almost insensitive to primordial non-Gaussianity, so that 
shell-crossing is not amplified and exact results can be ob- 
tained provided one uses the correct linear density thresh- 
old, rather than the usual one. We also propose a simple 
recipe to match the dependence on /nl of the large-mass 
tail while keeping the mass function normalized to unity. 
Then, in section [H we consider the two-point correlation 
of dark matter halos in real space, following the spirit of 
Kaiser (1984). Next, taking a Fourier transform we obtain 
the halo bias in Fourier space. Here, our aim is to show that 
one does not need to introduce free parameters to match 
the results of numerical simulations. Moreover, the nonlin- 
ear real-space expression is of interest by itself and it also 
allows one to check whether the "linearized" bias is valid 
on the range of interest. Finally, we conclude in section [5l 



2. Non-Gaussian initial conditions 

We focus in this paper on non-Gaussianities of the local 
type, where Bardeen's potential <& is of the form (|1]), with 
a Gaussian random field. On scales smaller than the Hubble 
radius, $ equals minus the Newtonian gravitational poten- 
tial and the Poisson equation gives in Fourier space (Slosar 
et al. 2008) 



^i,(k, z) = a(fc, z)$(k) with a{k,z) 



2c^eT{k)D{z) 

Sr^niO^O 



,(3) 



where 5l is the linear matter density contrast, T{k) is the 
transfer function and D[z) is the linear growth factor, nor- 
malized as D{z) (1 + z)^^ at high redshift. Unless stated 
otherwise, we normalize the Fourier transform as 



(5l(x) = / dke"' ''5L(k). 



Note that we define /nl by applying Eq.([T]) at early times 
(i.e. z = cx)), which is sometimes called the "CMB conven- 
tion" , whereas some authors first linearly extrapolate <i> at 
z = ("LSS convention"). Thus, both conventions are re- 
lated by f^^^ = -0(0)/^^ (PiUepich et al. 2009). Then, 
defining the time-dependent Gaussian field x by 

X(k, z) = a(A:,z)0(k), (5) 
we can write the linear density field at redshift z as 



^L(k) = X(k)+ / dkidk2Jl3(ki+k2-k)/^L(ki,k2) 



X X(ki)x(k2 



with 



local type : /NL(ki, k2) = /i 



NL 



Q:(ki + k2) 



(6) 



(7) 



a(fci) a{k2) ' 
This reads in real space as 

'5l(x)=x(x)+ j dxidx2/NL(x;xi,X2)x(xi)x(x2), (8) 
where the kernel /^l(x;xi,X2) only depends on the two 



vectors {xi — x, X2 — x}, 
/nl(x;xi,X2) = /nl(xi -x,X2 -x). 



(9) 



as long as the system remains statistically homogeneous, 
as for the local model ([T]). The real-space and Fourier-space 
kernels are related by (note the different normalization from 

O) 

/^l(xi,X2) = j ^M^e"'--+*--/^Jki,k2). (10) 

The relationships (O and ([5]) describe any homogeneous 
model where the Bardeen potential can be expressed as the 
sum of linear and quadratic terms over some Gaussian field. 
Thus, our analytical results also apply to other "/NL-type" 
models than the "local" one (|T]). 

For /nl = we recover Gaussian initial conditions, 5^ — 
X, with a linear density power spectrum 



(X(ki)x(k2)) =<5z,(ki+k2)PL(fcl), 

and a two-point linear density correlation 
Cl(xi,X2) = (x(xi)x(x2)) 



(11) 



= 4tt dkk^PL{k) 



sin(fc|x2 — xi 
fc|x2 - Xi| 



(12) 



As usual, it is convenient to introduce the smoothed linear 
density contrast, Xg(x), within the sphere of radius q and 
volume V around position x. 



X?(x) 



dx' 



X(x- 



dke*''x(k)T^(fcg), (13) 



with a top-hat window that reads in Fourier space as 



(4) W{kq)^ l 



i(A;g) — kg cos{kq) 



(14) 
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Then, in the linear regime, the cross-correlation of the 
smoothed linear density contrasts at scales qi and 172 and 
positions xi and X2 = Xi -I- x reads as 

^9i,92(^) = (Xgi(xi)Xq2(xi +x)) 



An / dkk^PL{k)W{kqi)Wikq2 



sm{kx) 



(15) 



In particular, aq — (Jq q(0) is the usual rms linear density 
contrast at scale q. Then, in the non-Gaussian case we de- 
fine the initial conditions by the same power spectrum (fTTj) 
for the field x and we vary the parameter /nl (for the local 
type (O). We still define the variance as in Eq. p3|) from 
the Gaussian field x- 

In the following sections, where we define dark matter 
halos as spherical overdensities, we shall need the average 
of the kernel /^l(x; xi , X2) over spherical cells, weighted by 
the linear correlation (jl2p . Thus, omitting the superscript 
S for simplicity, we define the quantity 

/,,,„,(x;x„X2) = X^X^^X^^/dq'idqi 

X /^L(q; q'l, q2) c^L(q'i, qi) CL(q'2, q2), (i6) 

where the spheres of volumes V, V\ and V2, and radii q, q\ 
and ^2, are centered on the points x, Xi and X2. In Eq. (fT6|) 
the coordinates q'j^ and qj are integrated over all space. In 
terms of the Fourier kernel this reads as 



/9;«i,g2(x;xi,X2) = y dkidk2/NL(ki,k2)PL(fci)FL(fc2) 

X W{k^qx) W{k2q2) W^(|ki + k2|q) 

X giki.(xi-x)+ik2.(x2-x) ^^7) 

Thanks to statistical homogeneity and isotropy, the kernel 
/^L(ki,k2) only depends on the lengths /ci, fc2, and on 
the angle between both vectors. Then, for spheres that are 
centered on the same point (i.e. x = xi = X2), Eq. (fT7|) 
simplifies as 

r-OO 

x = Xi=X2: /,;gi,,2 = SttM Ak^k\Pl^(k^)W(k^q^) 



X y Ak2kiPL{k2)W{k2q2) J dtiW{kq)f^^{ki,k2,^i), 

(18) 

which does not depend on the position x of the sphere. Here 
wc defined fi — (ki.k2)/(fciA;2), and k"^ = k\ + k2+ 2kik2^- 
On the other hand, when we consider the two-point corre- 
lation of dark matter halos, the three spheres in Eq. (fTB)l are 
chosen among two possible spheres Va and VJ,, separated by 
a distance x. Then, we need the two quantities, 

/•oc 

V = Va, Vi^V2=Vb: fa-bb{^) = / dfci fc?PL(fcl) 

Jo 

1*00 

XW(kiqb) / dk2klPL{k2)W{k2qb) 

Jo 

X / dMmfc9a)/^L(fcl,fc2,M)^^^, (19) 



kx 



and 



/■oc 

Vl =Va,V^V2^Vb: h;ab{x) = StT^ / dkl kjPLiki) 

Jo 



X^(fcl9a)^^^T^ / dk2p2PLik2)W{k2qb) 



kix Jo 
X / dnW{kqb)f^i^iki,k2,^i). 



(20) 



Their Fourier transforms with respect to the separation x 
read as 

fa-bbik) - / dkidk2(5l5(kl+k2-k)Pi(fci)Pi(fc2) 



xWikiqb)Wik2qb)Wikqa)fL0^lM), (21) 

and 

fb-abik) = PL{k)W{kqa) J dkl PL{ki)W{kiqb) 

xW^(|ki+k|<7;,)/^L(k,ki). (22) 

3. Mass function of dark matter halos 

We now extend the analysis of Valageas (2009b) to obtain 
the mass function of dark matter halos for non-Gaussian 
initial conditions. 

3.1. Rare-event saddle point for the density distribution 

In a fashion similar to Valageas (2009b), we note that the 
exponential falloffs of the large-mass tail of the halo mass 
function n{M), and of the overdensity tail of the linear den- 
sity contrast distribution Vhi^Lq), can be exactly obtained 
from the constrained maximum, 



rare events : Vl{Sl) 



max e 

{xl<i]\SL,=SL} 



-kx-c: 



(23) 



Here we introduced the probability distribution VLi^Lq) of 
the smoothed linear density contrast within the sphere of 
radius q, which we can take centered on the origin, 



d\<i~5L0^)W{kq). 



(24) 



In the Gaussian case, /nl = 0, we simply have 5^ = Xg, as 
defined in Eq. (fT3|l . and the probability distribution Vl is 
a Gaussian. Then, as stressed in Valageas (2009a,b), in the 
limit of rare events (e.g., at fixed density contrast S^q in the 
large-scale or large- mass limit q — > 00), the tails of the dis- 
tribution Vl are obtained from Eq. (|23p , where we maximize 

the statistical weight e"^-^''-'^ -^^/^ of the Gaussian field X: 
under the constraint that the smoothed linear density con- 
trast 5Lq is equal to the value of interest. Equation (^5)1 
only gives the leading-order exponential falloff. Subleading 
terms, such as power-law prefactors, may be obtained by 
expanding around the saddle-point Xi within a steepest- 
descent method. This also provides the probability distri- 
bution V{5) of the nonlinear density contrast 5 at scale r in 
the quasi-linear regime, as well as the cumulant generating 
function (Valageas 2002a, b; 2009a, bfl. As pointed out in 



^ As shown in Valageas (2009a), for the closely related adhe- 
sion model, where the same procedure can be applied, one can 
explicitly check that the asymptotic results obtained by this ap- 
proach agree with the complete distribution V{5) that is exactly 
known for two cases (Brownian and white-noise linear velocity 
in 1-D, corresponding to a power-law linear density power spec- 
trum with n = —2 and n — 0). 
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Valageas (2009b), this also gives the large-mass tail of the 
halo mass function, if halos are defined as spherical over- 
densities with a fixed nonlinear density threshold Sr- Then, 
the halo mass M and radius r are related to the Lagrangian 
radius q and linear density threshold Slq through 



the parameter /nl)- At order zero we recover the Gaussian 
saddle-point, 



(l + (5,)r^ with Sr = T{SLq), 



and 



M 



47r 



-9 , 



(25) 



(26) 



A(o) 



Sl 



(32) 
(33) 



At first order we obtain 
X 



where the function T describes the spherical collapse dy- 
namics. This holds as long as shell-crossing has not ex- 
tended beyond radius r, so that the usual spherical col- ^2 
lapse dynamics at constant mass is vahd. This yields an A^^^ = "3 -| fq;qq, 
upper bound 6+ for the nonlinear density threshold that '^q 
can be used to define halos to take advantage of the exact 
asymptotic tail (I23p . For the preferred ACDM model this 
gives 6+ - 200 for M - lO^'^h'^MQ up to S+ ~ 600 for 
M ~ lO^^h^^AlQ (the dependence on mass is due to the 
change of slope of the matter power spectrum Pi(A;)). 
In the Gaussian case, the constrained weight (j23p simply 

reads as e"*^'^^^'^'-', and we recover the exponential tail of 
the usual Press-Schechter mass function (Press & Schechter 
1974), except that the standard threshold 5c — 1.686 must 
be replaced by Sl = :F-\5), with 5l ^ 1.59 for S = 200. In 
the non-Gaussian case, that is for /nl 0, we can compute 
the weight l|23p by using a Lagrange multiplier A. Thus, we 
define the action iS[x, A], 



<^^(q)-A(^)/^^C.(q,q') + 2A(°)^/^^ 

X y"dqidq2/NL(q';qi,q2)C'L(q,q2)CL(qi,q'i), (34) 



(35) 



where the quantity fq.qq is given by Eq.([TH]) 



(27) 



where <5Lq[x] is the nonlinear functional that affects to the 
initial condition defined by the Gaussian field x the linear 
density contrast 5Lq within the sphere of radius q, obtained 
through Eq.([8|). Then, we must look for the saddle-point of 
the action (|?7|) . with respect to both x and A. Thus, dif- 
ferentiating the action with respect to x(q) a-nd mul- 
tiplying by the operator Cl gives 




X(q) = ACL(q,q') 



V6 



Lq 



-Dxici' 
whence, using Eq.®, 
dq' 



(28) 



Fig. 1. The radial profile ((36|) - ((37)l of the linear density 
contrast SLq' of the saddle point of the action cS[x, A]. We 
show the profiles obtained with a ACDM cosmology for the 
masses M — 10^^ and 10^^h~^MQ. A larger mass corre- 
sponds to a lower ratio 6Lq' I^Lq at large radu (f jq > 1. We 
show our results for the Gaussian case (solid line) , positive 
/nl (dashed line) and negative /nl (dotted line), for the 
local model ([7]). 



x(q) 



A /^-^CL(q,q') 



2A 



dq' 



dqidq2 



From Equations ([52|) - ([35)l we also obtain the radial lin- 
ear density profile of the saddle-point, up to first order. 



X /NL(q';qi>q2)C'L(q,q2)x(qi)- 



(29) .(0) 



Differentiating the action 1^7]) with respect to A gives the 
constraint 



= SLqix], 
whence 

dq 

-1/ 



(30) 



'Lq' 



5(1) 



^2 

~2 q^q' 



fq'\qq + '^fq.q 



q-q' 



fq;( 



(36) 



(37) 



x(q) + ydqidqa /NL(q; qi, q2) x(qi) x(q2) 



(31) 

Next, we solve the system (l29l) - ((3T]) as a perturbative series 
over the non-Gaussianity kernel /^^ (i-S- over powers of 



We can check that at q' = q it verifies the constraint pO|) . 
and at order zero we recover the Gaussian profile (Valageas 
2009b). Equations ([5S|) - ([57)l give the integrated density pro- 
file, that is, Sj^q' is the mean linear density contrast within 
the Lagrangian radius q' . The local linear density contrast 
at radius g', 5L{q'), is given by 



iq'^Liq') = ^{q'^hq'), 



(38) 
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whence 



3.2. Mass function 



SI 



(39) 



(40) 



where fo-qq{q') and fq-oq{q') are given by Eas. (fTg)) - ((^ . 

We show in Fig. [T] the integrated linear density pro- 
file ([55)l - ([57|) obtained for the masses M = 10^^ and 
IO^^H^^Mq, for a ACDM cosmology. The dependence on 
mass is due to the change of slope of the linear power spec- 
trum with scale. We plot our results for the Gaussian case 
(/nl = 0, sohd hues), large positive /nl (/nl = 2 x 10^ 
and /nl = 2 x 10**, dashed lines) and large negative /nl 
(/nl = — 2 X 10^ and /nl = — 2 x 10**, dotted lines), for the 
local model ([7]). Thus, we can see that a positive /nl in- 
creases the relative density (i.e. with respect to the density 
at radius q) both at small and large radii. The very large 
values of /nl required to be able to distinguish the curves 
in the figure imply that for realistic cases (|/nl| < 100) the 
perturbation of the density profile is very small. Therefore, 
the values of the upper boundary 5+ , which marks the on- 
set of shell-crossing, obtained in Valageas (2009b) for the 
Gaussian case remain valid up to a very good accuracy. 
We can note that to obtain a similar deviation from the 
Gaussian profile we need a larger value of the parameter /nl 
at smaller scale. This can be understood from the expres- 
sion ([7]), which scales as ~ fNh/a{k) oc f^i,/{k'^T{k)) 
and grows as fc~^ at very large scales. The same behavior 
(i.e. a higher sensitivity to local-type non-Gaussianity at 
large scales) is obtained for the bias of dark matter halos, 
see Eq.([2]) above and section 2] below. 

Next, we define the constrained weight ([^5]) as 



(41) 



at the relevant saddle-point x- This gives, up to first order, 



. .'-X--^, (42) 

r« = x(°).c,~^x« = (43) 

Therefore, the tails of the probability distribution ([25)1 read 
as 



I 4 



(44) 



where we introduced the skewness of the linear density con- 
trast, at first order over /nl, 



{^Iq? 



(1) 



= 6 



fq;qq 



(45) 



Indeed, from Eq.® we have at first order 



3 \(1) 



Lql " 1^ y3 



dqidq2dq3 



I dqldq'2/^L(q3;q'iq2) 



x(x(qi)x(q2)x(qi)x(q2)) 

— ^ fq;qq- 

Hereafter, since S^^^ = 0, we simply note S3 = S^^^\ 



(46) 
(47) 



Following the Press-Schechter approach (Press & Schechter 
1974), the mass function that is obtained in the Gaussian 
case from the probability distribution 'Pi (Jig) reads as 



n(A/,/NL = 0)dA/=^/(z.)^, 



(48) 



with (using the subscript "PS" to distinguish the Press- 
Schechter prediction) 



/ps(i^) 
and 



a{M) ' 



(49) 



(50) 



Here a{M) = cr^, where the Lagrangian scale q is related to 
M by Eq.lHSl). As stressed in Valageas (2009b), the linear 
threshold 5^ in Eq. ((50|) must be defined as !F~^{5), as in 
Eq.(Il5l), which gives 5l ^ 1.59 for 5 = 200. Then, the 
exponential tail of (|49p is exact, 

v^^x,: ln[/(i.)] ^ (51) 

but the power-law prefactor and the low-mass tail of (|49|) 
have no reason to be valid (and numerical simulations in- 
deed show that they are not exact). Then, in order to match 
numerical simulations, one needs to use fitting formulae. 
One such fit to simulations, which obeys the exact tail ([5T|) . 
is (Valageas 2009b) 

f{v) = 0.5 [(0.6 vf-'' + (0.62 y f-''] e-"' . (52) 



Both mass functions ([49]) and (|52|) satisfy the normalization 



1, 



(53) 



which ensures that all the mass is contained in such halos: 
Mn{M)dM ^-p^. (54) 



In the non-Gaussian case, that is /nl 7^ 0, we may es- 
timate the halo mass function by multiplying the Gaussian 
one by the corrective factor obtained in Eq. ((44|l . 



/nl ^ : n(Af,/NL) = n{M,Q) e^'''"-/'^^<\ 



(55) 



As explained above, this yields the exact large-mass tail (up 
to first order over /nl) but it is not expected to hold for 
the low-mass tail. In particular, this gives a non-Gaussian 
mass function that does not obey the normalization (|54p . In 
order to satisfy Eq. ((54)) . while keeping the large-mass tail 
of Eq. ([55)1 , a simple procedure is to make use of the scaling 
and of the normalization ([55]) . Thus, modifying the 
relationship ([SD]) as (see Eq. (|ii)) ) 



1 - ^53(M), 



a{M) 

we may use for the non-Gaussian mass function 
/nl ^ : n(Af, /NL)dM = ^/(m)^, 



(56) 



(57) 
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where we use the same scaling function / as for the 
Gaussian case This recovers the large-mass tail 

and satisfies the normalization ((54)l . Note that Afshordi 
& ToUey (2009) have recently proposed a similar rescal- 
ing as ([55)) . which they re-interprate as a modified effective 
variance u^. However, to obtain such a rescaling they use 
some approximations, such as the decoupling between the 
local values of the fields 0(x) and in the spirit of a 

peak-background split approximation, whereas the deriva- 
tion presented here is asymptotically exact. Equation (|57p 
yields for the ratio of both mass functions at fixed mass M, 



n(M, /nl) ^ f{^J■)/fJ■ d^i 
n(M,0) ^ f{iy)/iy dv' 

Using Eas. ((50l) . ((56)) . we obtain 



d^; = V'-y^^ 



6a/1 



dS^ 
da ■ 



(58) 



(59) 



If we use the Press-Schechter mass function P^ . Eqs. ([57|) - 
(|59|) give back the result obtained in Matarrese et al. (2000), 
except for the fact that we use 5l — 1.59 instead of 1.686 
as explained above. 

Note that Matarrese et al. (2000) also use a saddle-point 
approach to derive the halo mass function for non- Gaussian 
initial conditions, expanding at linear order over /nl (or 
5*3), so that their computation is equivalent to the one de- 
scribed above. However, since they first compute the cumu- 
lant generating function {W{X) in their notations, or Lp{y) 
in Valageas 2009b) , the constraint ((30)) is expressed through 
a Dirac function, written as an exponential by introduc- 
ing an auxiliary variable A, so that they can first integrate 
over the Gaussian field x s-nd next expand over /nl the ex- 
pression obtained for W{X). The somewhat simpler method 
described in this article has the advantage of simultane- 
ously giving the density profile ((55)) - ((57)) of the underlying 
saddle-point. As explained above in Fig.)!] this allows us to 
check that realistic amounts of primordial non-Gaussianity 
have a negligible effect on this profile, so that the onset of 
shell-crossing appears for almost the same nonlinear density 
threshold This ensures that the rare-event and large- 
mass tails ((31]) and ((55)) are exact (at leading order), as 
long as halos are defined by a nonlinear threshold 5 ~ 200 
below the upper bound i5+. The simplicity of the method 
presented in this article also allows a straightforward ap- 
plication to two-point distributions, as shown in section (D 
below, or to more complex primordial non-Gaussianities, 
that may involve several fields or higher-order polynomials 
as in Eq. ((95)) below. 

Another approach presented in Lo Verde et al. (2008) 
is to use the Edgeworth expansion, that writes the proba- 
bility distribution Vl{5l) as a series over the cumulants of 
the non-Gaussian variable 8l- In practice, one truncates at 
the lowest order beyond the Gaussian, that is, at the third 
cumulant {5\)c described by S'3. Thus, expanding the expo- 
nentials ((44)) and ((55)) one recovers the results of Lo Verde et 
al. (2008) at large v (i.e. low Cg). However, in the rare-event 
limit, where computations rest on firm grounds as explained 
above, the Edgeworth expansion does not fare very well. In 
particular, although we only derived the saddle-point x ^nd 
the argument P of the exponential up to linear order over 
/nl, see Eas. ((42t - ((43t . it is best to keep the exponential as 
in Eas. ((55)) or ((57)) . Indeed, in the rare-event and small- /nl 



limits, the tail ((55)) can be a good approximation even when 
the corrective factor Sz5\/ {Qa^) is much larger than unity 
(i.e. we can have the hierarchy P^") > P^^^ > 1 > P(^)). 

Since the low-mass tail ()49p does not match numerical 
simulations, it is sometimes proposed to keep the ratio ()58p 
given by the Press-Schechter mass function, and to multiply 
the fit from simulations of the Gaussian mass function by 
this factor (Grossi et al. 2007, 2009; Lo Verde et al. 2008). 
However, this procedure clearly violates the normalization 
condition ([53]). Therefore, we suggest to use ((55)) with the 
fitting formula obtained from Gaussian simulations, that is, 
to use Eq. ((57p . which automatically satisfies the normaliza- 
tion ()54p . However, there is no reason to expect that the 
low-mass tail can be exactly recovered by any such proce- 
dure, even though by construction it gives the right behav- 
ior for the Gaussian case, as the low-mass slope might also 
depend on /nl in some specific manner. 
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Fig. 2. The ratio n(M, /NL)/n-(Af, 0), of the mass func- 
tions obtained for /nl = ilOO over the mass function ob- 
tained for Gaussian initial conditions, as a function of M for 
several redshifts. The dot-dashed line labeled "e"^^" is the 



multiphcative factor ((55)) . the dashed hne labeled "PS" is 
Eq. ((58| with the Press-Schechter mass function ((44)) . which 
also corresponds to the result of Matarrese et al. (2000) (but 
with Sl ^ 1.59), the sohd line labeled "f" is Eq.((5H]) with 
the fitting function ((5^ . The data points are the results 
from the numerical simulations of Grossi et al. (2009). 

We compare in Figs. [2l [S] our results for the halo 
mass function with numerical simulations from Grossi et 
al. (2009). Note that they use the LSS convention for /nl, 
so that in terms of the CMB convention used in this article 
this corresponds to /nl — ±75 (Fig. ^ and /nl — ±151 
(Fig. [3]). We show the ratio of the non-Gaussian mass func- 
tion to the Gaussian one, as given by the multiplicative 
factor ()55p (dot-dashed line), the ratio ((55)) computed with 
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Fig. 3. Same as Fig. [2] but with /^p = ±200. 

the Press-Schechter mass function ([^^ (dashed line) or 
with the fitting function (|52p (solid line). In agreement 
with Grossi et al. (2009), we find that the ratio ([55)1 com- 
puted with the Press-Schechter mass function, which also 
corresponds to the result of Matarrese et al. (2000) (but 
with 5l — 1-59 instead of 1.686 as explained above), agrees 
reasonably well with simulations. Using the fitting func- 
tion ((52l) . or the simple multiplicative factor ((55)) . yields 
close results in this regime and also agrees with simula- 
tions. However, using Eg. ([55]) with the fitting function 
appears to give a somewhat better agreement with simula- 
tions, especially at low masses. This could be expected from 
the fact that this procedure ensures that the mass function 
is properly normalized (by contrast, the simple multiplica- 
tive factor (|55p is larger than unity and does not reproduce 
the crossing of both mass functions for ~ 1). 

Next, we compare in Fig. |3] our results with numerical 
simulations from Dalai et al. (2008). We again show the 
ratio n{M, fi^i)/n{M,Q) as a function of mass, but for a 
larger primordial non-Gaussianity, /nl — ±500. For /nl = 
500 (upper panel) we note as in Dalai et al. (2008) that the 
prediction of ([55)1 computed with the Press-Schechter mass 
function (|44|) (i.e. the result of Matarrese et al. 2000 but 
with 5l — 1.59) tends to overestimate the deviations from 
the Gaussian case (and the simple multiplicative factor (j55p 
fares somewhat worse). However, using Eq. ljSS]) with the 
correct Gaussian mass function ([5^ decreases somewhat 
this ratio (as in Figs. [H [3]) and provides a good agree- 
ment with the simulations. For /nl = —500 (lower panel) 
the match is not so good. However, in the latter case the 
agreement might improve at higher masses, where there 
are no data points but where the predictions fSS]) or ([58]) 
are asymptotically exact. Moreover, we can note that the 
discrepancy between theoretical predictions and numerical 
simulations is of the same order as the deviation between 
the different theoretical curves. Since the latter share the 




10" 101= 




M [17 -'Mq] 

Fig. 4. Same as Fig. [21 but with /nl — 500 (upper panel) 
and /nl = —500 (lower panel). The data points are the 
results from the numerical simulations of Dalai et al. (2008). 

same exact large-mass behavior (at the leading order given 
by the exponential cutoff (HH)), these deviations show the 
sensitivity of the mass functions to the details of the theo- 
retical prescriptions. Note that these power-law prefactors 
have not been rigorously derived (and are expected not to 
be exact, for all formulae used here). Therefore, the devia- 
tion between these theoretical predictions gives an estimate 
of the theoretical uncertainty for the ratio of the halo mass 
functions. Then, taking into account this theoretical uncer- 
tainty, we can see that the agreement with the numerical 
results is still reasonable. For practical purposes, when one 
tries to derive constraints on cosmology from observations 
of halo mass functions, it would be useful to consider several 
theoretical prescriptions in addition to the best prediction 
([55)) - ([5^ . as in Figs.[2]-Fig.[ll so as to take into account the 
theoretical uncertainty in the analysis. 

4. Bias of dark matter halos 

4.1. Two-cells saddle-point and real-space bias 

We now consider the bias of dark matter halos, or more pre- 
cisely their two-point correlation function. As in Valageas 
(2009b), following Kaiser (1984), since we identify rare mas- 
sive halos with positive density fluctuations in the linear 
density field, we first consider the bivariate probability dis- 
tribution, 'Pl{5li,5l2)^ of the linear density contrasts 5li 
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and 5l2 within two spheres Vi and V2, of radii qi and q2, 
and separated by the distance s. Note that we distinguish 
the Lagrangian distance s between the halos, measured in 
the linear density field, from their Eulerian distance x, mea- 
sured in the nonlinear density field. Indeed, since halos have 
moved through their mutual gravitational attraction these 
two distances are usually different. Proceeding as in sec- 
tion [23 to obtain the rare-event tails. 



Vl{5li,5l2) max e ^x-C^ 

{x[q]|<5-Lgl=l5l-l,l5Lg2=l5i2} 

we are led to introduce the action 

^[x, Ai,A2] = Ai - (5l,i[x]) + A2 ((5l2 - ^Lq2[x]) 
1 

+ 2 



(60) 



oX-t^Z'-X (61) 



which now involves the two Lagrange multipliers Ai and 
A2 . We again obtain the minimizer of the Gaussian weight 
ri2 — (x-C'^^-x)/2 of Eg. ([SO]) by differentiating the action 
<S with respect to Xj Ai and A2, and we solve these equations 
as a perturbative series over /nl- At order zero we recover 
the Gaussian terms 

x(")(q) = Ar/ ^C,(q,q,) + A^"Y ^C.(q,q2) (62) 

(63) 
(64) 



A 



(0) 



A, 



(0) 



crj Sli - 


C^?2 Sl2 


erf 


4 

-0-12 


crl 5l2 - 


cr?2 Sli 




-f^f2 


note af 





Eq. (fT5|) . At first order we obtain 
X«(q) = A«/ ^C.(q,qi) + A« / ^C.(q,q2) 

JVi ^1 JV2 ^2 

+2A("y^^ idq'idq^ /^Jqi; q;, q^) Ci(q,ql) X^°Hq2) 

+2A("y^ ^ /^L(q2; q'l, q2) ^^(q, ql) x'"^ (q2) 

(65) 

X['\alal - af^) = X^^^' [aU^fuu + f2;n) - Sa^h.n] 

+ 2Xf>h^°^ [^?2(2/2;12 + /l;22) - ^2'(2/l;12 + /2;ll)] 
+ Af ' [3a22/2;22 - a|(2/2;12 + /l;22)] , (66) 

and 

Xi'\alal - af^) = X^"^' [3alj,,n - <j!{2f,,u + /2;ii)] 

+ 2A(")a(") [^T?2(2/i;12 + /2;ll) - <jl{2f2a2 + /l;22)] 

(67) 



Af' [af2(2/2;12+/l;22)-3af/2;22] 



where we note for instance /i-22 — fa-bb{s) and /212 — 
h;ab{s), with Va = , Vb = V2, from Y^qs.^-^. Ncxt, 
the Gaussian weight ri2 = (x.C'^^.x)/2 reads at order zero 



p(o) _ Sli + ol - 2a\^ 5liSl2 
^ 12 — 



and at first order, 

= A« A^ + .f2Af ) + A« {al^X^^^ + .fA^")) 

+2Ar^^al+2ArAf (2/ia2 + /2;n) 

+ 2A("^A(°^'(2/2a2 + /l;22) + 2A(°)V2;22. (69) 

Therefore, in the rare-event limit the tail of the bivariate 
distribution (pD|) reads as 

-Ar 



Vl{6li,Sl2) - PLiSLi)rLiSL2)e- 
with 

Ar-ri2-ri -r2, 



(70) 
(71) 

where ri2 is given by Eas. ([55)l - (|M|) . and Fi and r2 by 
Eqs.(l42])-([43l) for each sphere Vi and V2. Next, following 
Valageas (2009b), we write the halo two-point correlation 
as 

l+ah,M2ix) - il + SLM{s))e-^^, (72) 

where the factor (1 + Slm{s)) models the effects associ- 
ated with the mapping from Lagrangian to Eulerian space. 
This is the local linear density contrast at radius s from 
a halo of mass M — max(Mi, M2) (to keep the symmetry 
Ml ^ M2), as given by Eqs.(I2ni)-(Iini). A sufficiently accu- 
rate approximation would be to use only the zeroth-order 
term (P^ . as shown by Fig. [2 but taking into account the 
correction (j40p brings no further difficulty. Here we approx- 
imated the nonlinear density contrast 6m by the linear den- 
sity contrast 5lm, since at large distance where 6m <C 1 we 
have 6m — Slm- Next, we must express the Lagrangian 
separation s in terms of the Eulerian distance x. Following 
Valageas (2009b), at lowest order where we consider each 
halo as a test particle that falls into the potential well built 
by the other halo, we obtain s as the solution of the implicit 
equation 



1 



As) 



>LM: 



(73) 



where 6LMi{s) is the linear density contrast within radius 
s of the halo of mass Mi, given by Eas. p6| - p7|) . At large 
separation, this relation can be inverted as 



A , SlMi jx) 6lm2 jx) \ 



(74) 



which provides an explicit expression for s. 

Finally, we define the real-space halo bias as the ratio 
of the halo and matter two-point correlations, 



^Mi.M2 (^) 



^Mi,M2(a^) 



(75) 



Since at large distance the matter correlation is within the 
linear regime, ^{x) ~ CTq ol-^); ^® ^-^^o write in this limit 



(76) 



which fully determines the halo bias from Eg. ((72|) . 

For equal-mass halos of radius q, defined by the same 
threshold 6l, the two Lagrange multipliers are equal, Ai = 
A2 = A, and Eqs.(l63l)-(l64l), simplify as 



2{alai-af2) 



(68) A(°) = 



6l 



a'' + al2 



(77) 
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and 



-3 Si 

(/2;11 + 2/l;12 + /l;ll), 



(78) 



where we note = <^q- This yields for the Gaussian weight 



12j 



.(0) 
12 



and 



SI 



t2 I ^2 ' 



'12 



(1) 



^2 SI 



(/2;11 + 2/l;12 + /l;ll) 



(79) 



(80) 



Then, the difference AT defined in Eg. ([TTjl . that measures 
the correlation between rare events in the linear density 
field, writes as 



and 



Ar(i) : 



(81) 



-2 SI 



(a2 



'12^ 



/2;11 + 2/l;12 + 



/1;11 



(82) 



We can check that AF vanishes for s — > oo, since all mixed 
quantities, ct^jj /i;12 and /2;ii, go to zero. 

As stressed in Politzer & Wise (1984), real-space for- 
mulae such as ((72|) . which are obtained in the rare-event 
{Sl/<^ ^ 1) and large separation {ai2 <C a) limits, do not 
assume that the exponent AF is small. In fact, as shown in 
Valageas (2009b), at large redshift one can probe a regime 
where this exponent is large, so that one needs to keep the 
nonlinear form ([7^ . which yields a nonlinear bias. As seen 
from Eq. (j81[) . this regime corresponds to very massive ha- 
los, ^l/ct — > oo, at fixed (small) ratio OYija. Nevertheless, 
in the regime where AF is small (i.e. in the large separation 
limit X — * oo), which covers the cases of interest encoun- 
tered at low redshift, we can expand the exponential in 
Eq. (|7^ . Then, at lowest order over the terms that vanish 
in the large separation limit, we obtain for equal-mass halos 



X ^ oo : ^M(a;) ~ (5lm(s) - AF, 
whence 



(83) 



si ' 



/o;,,(s) + 2/,;o,(s)-3%^a,%(s) 



/2;ll(s) + 2/ia2(s)-3^a,%(5) 



(84) 



Hereafter we call the bias 6^/(x) obtained from Eg. ([84]) and 
Eq. (|75|) as the "linearized" bias. 

We can note that, following the approach of Kaiser 
(1984), Matarrese & Verde (2008) also computed the effect 
of local-type non-Gaussianity H]) on the halo two-point cor- 
relation. Drawing on earlier work by Matarrese et al. (1986), 
who compute the n-point halo correlations by expanding 
the relevant path-integrals and next resumming the series 
within a large-distance and rare-event approximation, they 
obtain expressions of the form ((72)) . without the prefactor 
(1 -I- Slm), and ([M)) . without the terms in the first line. 



that arise from this prefactor, and the last term /i;ii in 
the second bracket. This latter term arises from the factor 
<j\i in the denominator of expression (j80p . associated with 
the effect of primordial non-Gaussianity on the one-point 
distribution l|43p . It can be seen as a renormalization of the 
Gaussian term (<5|,/o'^)ct^ ,j(s), since it shows the same scale 
dependence through the function (j\ q{s). The presence of 
such a term was already noticed in Slosar et al. (2008) in 
Fourier space, and Desjacques et al. (2009) pointed out that 
it needs to be included to obtain a good agreement with nu- 
merical simulations. This will give rise to the last term in 
Eq.([5n]) and the second term in Eqs.dH]) and ([Ml) below. 
On the other hand, the terms in the first line of Eq. (|M|) . 
associated with the prefactor (1 -I- Slm) in Eg. ((72)) . are due 
to the mapping from Lagrangian to Eulerian space, and 
the first bracket in Eg. ([M)) expresses the (small) effect of 
primordial non-Gaussianity on this mapping. Another dif- 
ference between Egs. ([7^ and and previous works is 
that (within some approximation) we take into account the 
difference between Lagrangian and Eulerian distances s and 
X. This can play a non-negligible role as seen in Fig. [5] be- 
low. 



X — 50 ii -JMpc 
z = 




1/ct(M) 

Fig. 5. The halo bias bM{x), as a function of cr(Af), at fixed 
redshift z — and distance x — 50/i~^Mpc. The solid lines 
"6" are the nonlinear theoretical prediction of Egs. ([72)) . ((73| 
and ([TG]) . for /nl = ±200 and /nl = (i.e. Gaussian case, 
intermediate line), while the dot-dashed lines "6^" are the 
linearized bias of Eq. ([84| . The upper dashed line "s = x" 
shows the result obtained in the Gaussian case by setting 
s = a; in Eg. (j72p . The points are the fits to Gaussian numer- 
ical simulations, from Sheth, Mo & Tormen (2001) (crosses) 
and Pillepich (2009) (circles). 



We show in Fig. [5] the halo bias bM{x), as a function 
of (t{M) at fixed distance x = 5Qh~^ Mpc and redshift 
z = 0, for /nl = ±200 and /nl = 0. We display both the 
nonlinear result (solid line) of Eq. (|72p and the linear result 
(dot-dashed line) of Eq. . In addition, for the Gaussian 
case we also display the bias obtained by setting s = x 
in Eg. ([72]) . As in Valageas (2009b), for the Gaussian case 
(/nl — 0) we obtain a good agreement with the fits to 
numerical simulations of Sheth, Mo & Tormen (2001) and 
Pillepich (2009). Moreover, we can see that it is important 
to take into account the displacement of the halos through 
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Fig. 6. The halo bias huix) as a function of distance x, at 
redshifts z = (upper panel) and z = 1 (lower panel), for 
several masses. We show the cases /nl = (solid hues), 
/nl = 100 (dashed hues) and /nl = —100 (dotted hues). 
The divergences at x ~ 120ft.~^Mpc are due to the fact that 
the halo and matter correlations do not change sign at the 
same distance. 



Fig. 7. The real-space ratio [A&M(a;, /NL)/&M(a;, 0)]/a:|o of 
the correction A6M(a;, /nl) = &M(a:, /nl) - bMix,0) to 
the Gaussian bias 6m (a;, 0), divided by the factor x^q with 
2^50 = a;/(50ft.~^Mpc), as a function of distance x. We show 
the cases /nl — 100 (upper hues) and /nl = —100 (lower 
lines) for several masses at redshifts 2 = (upper panel) 
and z = 1 (lower panel), from Ea. (|72p . Larger masses have 
a steeper slope at low x. 



Ea. (|73p . as the approximation 5 = 2: significantly overesti- 
mates the bias. We checked that using the simpler Eq. ([7i|) 
gives a close result to the one obtained with Eg. ([75)1 and 
also agrees with the simulations. Thus, for practical pur- 
poses it is sufficient to use Eg. ([74l) . We can see that for 
all cases shown in Fig. [5] the linear bias from Eq. gives 
results that are very close to the fully nonlinear expression 
([72]) . This justifies the use of such linearized expressions 
in this regime. This will be especially useful in section 
where we consider the bias of dark matter halos in Fourier 
space. Indeed, it is easier to take the Fourier transform 
of Eq. ([5^ , which allows to recover the results obtained in 
previous works. Thus, one interest of the real-space results 
l(7^ - ([M|) is to provide a check on whether linearized pre- 
dictions (i.e. where the halo correlation only involves the 
matter power spectrum at linear order) are valid. Then, in 
agreement with those studies (that were mostly performed 
in Fourier space and led to Eq.®, see Dalai et al. 2008; 
Slosar et al. 2008), and with Ea. ([84| . we can see that the 
deviation from the Gaussian bias, bpiix, /-ml) — bM{x,0), 



grows linearly with bM{x, 0) at large bias and has the same 
sign as /nl- 

Next, we display in Fig. [B] the dependence on the dis- 
tance X of the bias obtained for several masses at red- 
shifts z = and z = 1. More precisely, we show the ra- 
tio y/\^M{x)/^{x)\, since the halo and matter correlations 
do not change sign at the same point. We plot the cases 
/nl = ±100 as well as the Gaussian case /nl = 0. While 
the Gaussian bias is roughly constant over large scales, 
up to 100/i~^Mpc, in agreement with previous studies 
(e.g.. Mo & White 1996; Mo et al. 1997), the non-Gaussian 
bias shows a strong scale dependence, with a deviation 
from the Gaussian bias that roughly grows as x'^ up to 
^ 100/i~^Mpc. This agrees with the k~'^ behavior observed 
in Fourier space, see Eq.([2|) above (Dalai et al. 2008; Slosar 
et al. 2008) and Eq.dlSD below. 

In order to see more clearly the scaling of the real- 
space correction A6A/(a;, /nl) = &M(a;,/NL) - bMix,0) 
to the Gaussian bias, we show in Fig. [7] the ratio 
[AbMixjT<!L)/bMix,0)]/xlo, with X50 = x / (bOh-^Mpc) . 
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Fig. 8. The halo (/nl = 0, ±100) and matter (^) two-point 
correlations at redshifts z — (upper panel) and z — 1 
(lower panel). We show the curves obtained for the masses 
M = 1013 and 10"/i-iMn. 



We display the results obtained from Eg. ((72)) for several 
masses at z = (upper panel) and z — 1 (lower panel), 
for /nl = ±100- We can see that over the range 30 < a; < 
QO/i^^Mpc all curves roughly collapse onto one another. 
This means that over this range the real-space correction 
roughly scales as AbM{x) oc fNLbM{x,0)x^ , which roughly 
agrees with the Fourier-space scaling ([2]) (here we neglected 
any constant offset, such as the factor —1 in Eq.®, see the 
discussion of Ea. ([M)) below). In fact, it appears that our 
predictions scale more closely shown in Fig. [3 

than as a{l/x)~^, that would be suggested by Eqs.®, ([55]) 
(at these scales the transfer function already deviates from 
unity). This agrees with the behavior observed in Fourier 
space in Fig. [TT] below. Note that the masses shown in Fig. [7] 
span the range 1.2 < bM{x, 0) < 5.9 at x = 50/i~iMpc and 
z = 0, and 2.4 < bMix,0) < 15.4 at z = 1, so that the 
linear scaling with &M(a;,0) of the correction AbMix, f^i^) 
appears to be a good approximation. 

Below SO/i^^Mpc large masses show a steeper scale de- 
pendence for A6a/(x, /nl)- At very large distance, x > 
100/i~iMpc, the oscillations seen in Figs. [H] and [7| are due 
to the baryon acoustic oscillation. Indeed, since the halo 
correlation is not exactly proportional to cr§ p i^) i even in 
the Gaussian case and in the linear regime (for instance it 



involves the smoothing scale q, see Eg. ([5^ ). the baryon os- 
cillations seen in the halo and matter two-point correlations 
are not exactly proportional. This yields the non-nionotonic 
behavior seen in Figs. [Hand [7] around lOO/i^^Mpc. For the 
same reason, the halo and matter correlations do not ex- 
actly vanish at the same distance, which gives rise to the 
divergent spike at ~ 120h~^ Mpc. These features simply 
mean that at these scales it is no longer useful to work 
with the bias 6m, which only makes sense if the halo and 
matter correlations are roughly proportional. In this range, 
where the correlations show some oscillations and change 
sign it is no longer a good approximation to write the halo 
correlation in terms of the matter correlation multiplied by 
some slowly varying bias factor. Then, one rather needs to 
directly study the halo and matter correlations themselves. 

Thus, we compare in Fig. [5] the halo and matter two- 
point correlations. We focus on large scales to see how the 
baryon acoustic oscillation is modified when one uses mas- 
sive halos as a tracer of the initial matter power spectrum. 
In agreement with previous works (Desjacques 2008) we 
can see that the oscillation is strongly amplified for massive 
halos that have a large bias. This amplification still holds 
for significant primordial non-Gaussianity (/nl = ±100), 
although it appears to be slightly lower for positive /nl- 
Moreover, the peak of the oscillation shows no significant 
shift, so that a measure of its position appears to be a robust 
ruler to constrain cosmology, independently of the halo bias 
and of the primordial non-Gaussianity. By contrast, we can 
note that the distance at which the two-point correlation 
changes sign is not significantly modified as one goes from 
the matter to the halo correlation in the Gaussian case, 
but it is rather sensitive to the primordial non-Gaussianity. 
In particular, a positive /nl shifts this point to a larger 
distance. However, theoretical and observational error bars 
may be too large to use this effect to constrain /nl hi a 
competitive manner as compared with other probes. 

4.2. Fourier-space bias 

Rather than the real-space two-point correlation, recent 
works have mostly studied the effect of primordial non- 
Gaussianity on the halo power spectrum, where at lowest 
order the Poisson equation ^ directly gives an estimate of 
the form ^ for the deviation from the bias obtained with 
Gaussian initial conditions (Dalai et al. 2008; Slosar et al. 
2008). 

It is not convenient to take the Fourier transform of 
the nonlinear correlation ([7^ , but at moderate redshifts we 
have seen that the linearized form ([5^ provides a very good 
approximation. Then, if we also make the approximation 
s ~ X, which is valid at lowest order, the Fourier transform 
of Eq. readily gives the halo power spectrum as 



PM{k) ~ ^Pi^{k)W{kq) 



-3^PL{k)W{kq) 



SI 



/o;,,(fc)+2/,;0,(fc) 



^P^{k)W{kqf 



A 



/i; 



/2;ll(fc)+2/l;12(fc)-3^PL(fc)W^(fc(7) 



(85) 



where the quantities / are obtained from Eqs. pT|) - (l22|) . As 
discussed below Eg. ([M|) . the terms f2-ii{k) and 2/i;i2(fc) 
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in Eg. ([55]) were already obtained by Matarrese & Verde 
(2008), following also Kaiser (1984) by identifying mas- 
sive halos with rare fluctuations in the linear density field. 
Defining the Fourier-space bias as (note that this is not the 
Fourier transform of the real-space bias ([75)1 ) 



bUk) = 



PAijk) _ PMjk) 

P{k) - PL{k) - 



(86) 



where we used P{k) ~ Pbik) at low k for the matter power 
spectrum, we obtain h\j[k) from Eg. ([55]) . We can also ob- 
tain the bias of different-mass halos in a similar fashion, 
first expanding Eq. ((72|) and next taking the Fourier trans- 
form. In order to take into account the displacement of the 
halos (i.e. x ^ s), we can also make the approximation 



PAiik) 



a;\ 3 /a; , 



with X — 

k 



(87) 



where we use the explicit expression ([7^ for s{x). This 
expresses the fact that because of the displacement of mas- 
sive halos, which usually have come closer because of their 
mutual attraction (x < s), Lagrangian-space wavelengths 
~ s (i.e. measured in the linear density field) correspond 
to smaller Eulerian-space wavelengths ^ x (i.e. measured 
in the nonlinear density field). This also follows the spirit 
of the Hamilton et al. (1991) ansatz, translated in Fourier 
space in Peacock & Dodds (1996). 

We can note that for large negative /nl the halo power 
spectrum ([85)) can become negative. This is due to the fact 
that we looked for an expression of the halo correlation, 
or of the halo power spectrum, at linear order over /nl 
as in Eg. (|85p . However, by definition the power spectrum 
PAf (fc) must be positive. Then, in cases where the expres- 
sion (I85p turns negative one should consider higher-order 
terms over /nl, which would ensure that the power spec- 
trum remains positive. Nevertheless, since such high-order 
terms are beyond the scope of this article, wc consider below 
the following simple procedure that ensures that Puik), or 
the squared bias 6^(fc), remain positive. Up to linear order 
over /nl, the bias ([551) reads as 



^M(fc, /nl) 



/PM(fc,0) + APM(fc,/NL) 



&Af(fc,0) 



P{k) 

^ l APM(fc,/NL; 

2 PAf(fc,0) 



where APMikJ^h) = PAiikjNh) - PM(fc,0) is the devi- 
ation from the Gaussian halo power spectrum and we ex- 
panded the square-root. Then, we may use the last expres- 
sion ([88)) as the prediction for the halo bias. This amounts 
to make the transformation 



M = 2xl0iah-'Mg 
z = 0.2B 




k [hMpc-i] 




k [hMpc-i] 





k [iiMpc" 



k [JiMpc- 




k [AMpc-'l 



Fig. 9. The correction to the halo power spectrum 
due to primordial non-Gaussianity. We show the ratio 
PM{k, hiMPhiik, 0) - 1 for /nl = 100 (upper curves) and 
/nl —100 (lower curves). The solid and dashed curves 
that are almost indistinguishable are Egs. ([55)) and ([57)) . 
The dot-dashed curves, that are above the solid curves at 
low fc, correspond to Eg. ([89)) which ensures that the halo 
power spectrum is always positive. The theoretical predic- 
tions are for M = 2 x 10^^/i~^Mq and the data points are 
the numerical simulations of Desjacques et al. (2009), for 
Af > 2 X IQ^^h-^Mp,. 



+ 2^ (/2;ll(fc)+2/l,12(fc) 



(90) 



^\i{k^ /nl) 



6i,(fc,0) 



1 

^+2 



bli{k,h-L) 
6|,(fc,0) 



1 



(89) 



where the two sides only differ by higher-order terms over 
/nl (as 1 -)- e ~ (1 -I- e/2)^ up to hnear order over e) and 
the right side is always positive. 



If we take the limit of very rare events, that is cr„ 







in Eg. ([55)) . we can only keep the last two terms (note that 



/ oc 0-2 ggg Eq.dHD) below). 



: PM{k) ~ PL{k)W{kqf 



SI 



At low k, with W{kq) ~ 1, and at linear order over /nl, 
the square root of Eq. ([86)) gives with Eq. ([90)) the bias 



0, A: ^ : &M(fc) 



Sl 

rr2 



+ 



dl f2Mk) + ViMk) 
< PUk) 



(91) 



The first term, S^/a'^, is the result obtained by Kaiser 
(1984) for rare massive halos (i.e. with a large bias) for 
Gaussian initial conditions. It is interesting to note that 
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0.01 0.1 
k [JiMpc-'l 



Fig. 10. Same as Fig.El but for M ^ 1.5 x lO^/i^^ A/q. The 
data points are the numerical simulations of Desjacques et 
al. (2009), for lO^^ < M < 2 x W^h'^MQ. 



2 X VS^^h ^Mq, and we show in Fig. [TO] our results for 
M = 1.5 X VS^^hr^M^^ whereas the data points are for 
10" < M < 2 X IO^/i-^Mq. The predictions ^ (solid 
lines) and ([87|l (dashed lines) are almost indistinguishable 
in this regime and agree reasonably well with the simula- 
tions, except at low fc for /nl — —100 where they give a 
negative halo power spectrum. The dot-dashed line, which 
corresponds to Eg. ([55)1 . gives a much better fit to simu- 
lations at low k for negative /nl, as could be expected 
from the fact that it always gives a positive halo power 
spectrum. However, a priori one should not give too much 
weight to this improved accuracy in this regime. Indeed, as 
is clear from Eqs. ([55)l - (|5^ . the solid and dot-dashed curves 
in Figs. [MTUlonlv differ by terms of order and beyond. 
Since all our results have been derived at linear order over 
/nl one can expect that Eq. ([59]) does not include all terms 
of order /nl- Then, although for practical purposes it is 
better to use Eq. ([89l) in this regime (i.e. negative /nl at 
low fc) it is still useful to also consider Ea. ([55)) . as the de- 
viation between both predictions should give an estimate 
of the theoretical uncertainty. Let us point out here that, 
in contrast with some previous approaches (e.g., Grossi et 
al. 2009; Desjacques et al. 2009), the good agreement with 
numerical simulations shown in Figs. [9] and [10] is obtained 
from Eq. ([55)) without any fitting parameter (such as the 
rescaling parameter q in Grossi et al. (2009) or the mass 
function parameters in Desjacques et al. (2009)). 



the second term gives a scale-independent correction to the 
Gaussian bias. The presence of such a term was already no- 
ticed in Slosar et al. (2008) and Afshordi & ToUey (2009). 
As stressed in Desjacques et al. (2009), taking this term 
into account is required to obtain a good match to numer- 
ical simulations. We checked that this is indeed the case to 
match the numerical results in Figs. and 1101 below using 
Eq. (|85p . The last two terms depend on the wavenumber fc. 
There is an additional term, that we neglect in this paper, 
that arises from the dependence of the matter power spec- 
trum on /nl (i-e. the denominator in Eq. (|86p ). However, 
since this term is much smaller than the other ones we dis- 
regard it here (see Desjacques et al. 2009). For the local 
non-Gaussianity ([7]), we have from Eqs. iPH) - 



0: /2;ii(fc) - /NLa(fc) / dkiPi(fci) 



/l;12(fc) ^ /l 



NL ■ 



a{k) 



which yields from Eq. ([9T 



0, fc^O: 6M(fc) ~ 3|-3-|/i;ii+/i 



/NL 



261 



(92) 



.(93) 



Thus, we recover the fc~^ dependence at low fc brought 
by the local-type non-Gaussianity llj, through the l/a(fc) 
factor in the last term. The second term in Eq. ([93|) is the 
constant shift due to the non-Gaussianity noticed above. 

We compare in Figs. [9] and [10] our results for the 
Fourier-space bias h\j(k) with numerical simulations from 
Desjacques et al. (2009). Since the dependence on mass 
is rather weak, we show in Fig. [9] our results for M — 
2 X 10^^/i~^Mo, whereas the data points are for M > 




Fig. 11. The Fourier-space ratio 

A6M(fc, /NL)/Ao5M(fc,/NL) of the correction 
AbMik, /nl) = bM{k, /nl) — buik, 0) to the Gaussian bias 
6Af(fc, 0), divided by the scaling factor Aq^m defined in 
Eq. (|M)) . as a function of wavenumber fc. We show the case 
/nl = 100 for several masses, M = 10" (blue), lO" (red) 
and 10^^/i~^Mq (green), at redshifts z = (left panel) 
and z — \ (right panel). The curves labelled "s = x", 
"s > x" and "Pa/ > 0" correspond to Eqs.dSS]), dSTJ) and 
(IM)) respectively. 

As in Fig. [7] in order to see more clearly the scaling of 
the correction Ab^iik, /nl) = &M(fc, /nl) — i>M(k, 0) to the 
Fourier-space Gaussian bias, we show in Fig. [11] the ratio 
AbMik, /NL)/Ao&j\/(fc, /nl), where we now define 



Ao&Af(fc,/NL) = /NL&J\/(fc,0) 



2Sl 
a{k)' 



(94) 



The correction obtained from a simple peak-background 
split argument rather gives A6A/(fc, /nl) = /NL[&Af(fc,0) — 
l]{26L/a{k)), as in Eq.©, see Dalai et al. (2008); Slosar 
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et al. (2008). However, in our formalism both the leading 
term 5\/<7^ (third term in Eg. ([55]) ) and the subleading term 
5i/(T^ (first term in Eq. (j85p ) are modified by primordial 
non-Gaussianity (the terms in the brackets that follow). 
The term (5l/ct^ arises from the Lagrangian to Eulerian 
space mapping (the prefactor {1 + 5lm) in Eg. ([7^ ') and it 
corresponds to the factor 1 in the more usual Eq.(l2]), that is 
not modified in the simplest model. However, in general we 
can expect non-Gaussianities to also affect this Lagrangian 
to Eulerian space mapping, and our model gives an esti- 
mate of this effect through Eq. (|^0)) . Therefore, we scale 
Ao6M(fc, /nl) with hM{k, 0) rather than with [6m(^, 0) — 1] 
in Eg. ((94)) . We show in Fig. [7] the results obtained at red- 
shifts z = and z = 1 for several masses. The curves la- 
belled "s = x", "s > x" and "Pm > 0" correspond to 
Eqs.dHH), dHZl) and (HH) respectively. We can see that the 
collapse of the curves obtained for different masses is not 
exact, as could be expected since Eg. (|85p is more complex 
than Eg. ([94|) . but the scaling ([94)l still captures most of 
the dependence on the halo mass M . In fact, most of the 
dispersion seen in Fig. [7] does not arise from the different 
masses but from the various approximations ((85)) . ((87)) and 
(|5^ . This should provide an estimate of the theoretical 
uncertainty. As in Figs. [51 and [TUl Eg. (|59)) . which involves 
higher-order terms over /nl to ensure that Pmik) is always 
positive, rises above Eg. (|55)) at low k [k < O.Ol/i Mpc~^) 
and is indistinguishable at higher k, whereas Eq. (|57)) rises 
above ([85]) at high k {k > O.Olh Mpc~^), where the cor- 
rection A6M(fc, /nl) is quite small, and is indistinguish- 
able at lower k. Approximations ((M)) and (|57)) are likely to 
be most accurate at low and high k respectively, see also 
Figs. [H] and [Tni In any case, we can see that the scaling 
with wavenumber, Abiyiik, fjsih) ~ ^/ct{k), is slightly bro- 
ken, although this still provides a good approximation. In 
agreement with the numerical simulations of Desjacques et 
al. (2009), we can see that the ratio Ab/ AqB is suppressed 
at high k, which means that the correction A6j\/(fc, /nl) de- 
creases slightly faster than l/a{k) over the range 0.01— O.l/i 
Mpc""'^. In fact, in agreement with the real-space behavior 
seen in Fig. [7l it appears that over this range the correc- 
tion A6Af(fc, /nl) scales slightly more closely as than 
as l/a{k), although none of these two behaviors is exact. 



5. Conclusion 

We have shown in this article how to extend to non- 
Gaussian initial conditions the computation of the mass 
function and of the bias of dark matter halos presented in 
Valageas (2009b) for the Gaussian case. This relies on a 
saddle-point approach, that allows to derive the large-mass 
asymptotic tails of the quantities of interest from the statis- 
tical weight of the initial conditions, supplemented by ad- 
ditional nonlinear constraints. Then, focussing on the case 
of "/NL-type" primordial non-Gaussianity, where the linear 
gravitational potential can be written as the sum of lin- 
ear and quadratic terms over an auxiliary Gaussian field x, 
we explained how to obtain the relevant saddle-points as a 
perturbative series over the nonlinear parameter /nl- This 
method is very general and it applies to any case of small 
primordial non-Gaussianity, where Bardeen's potential (or 
equivalently, below the Hubble radius, the gravitational po- 



tential or the density field) can be written as a polynomial 
over a Gaussian field x 

k 

'^^X + T.fdx'-ix')), (95) 

where the nonlinear parameters fi are small. Then, follow- 
ing the method presented in sections [21 [H the one-cell and 
two-cell saddle-points (associated with the one-point and 
two-point density distributions, whence the halo mass func- 
tion and bias) can be computed as a perturbative series 
over the coefficients fi (which need not be of the same or- 
der). This includes in particular the case of a cubic term 
gNLX^- Of course, this also extends to the cases where the 
coefficients fi are not mere numbers but convolution ker- 
nels and to several Gaussian fields Xj (which can have a 
nonzero cross-correlation). In all such cases, the large-mass 
asymptotics are set by the statistical weight e~^^'''~^^^ ■Xj)/^ 
of these Gaussian fields, taken at the saddle-point associ- 
ated with the maximization of this weight under appropri- 
ate nonlinear constraints, that express the mapping from 
Xj to the relevant quantity, such as the nonlinear density 
within a spherical cell. 

Focussing on the case of local-type primordial non- 
Gaussianity, we described how to obtain, up to linear or- 
der over /nl, the one-cell saddle-point associated with the 
probability distributions Vl{Sl) and P{S) of the linear 
and nonlinear density contrasts within spherical cells. This 
gives the quasi-linear limit of these distributions as well as 
the large-mass exponential falloff of the halo mass func- 
tion. An advantage of our method is that it allows us to 
explicitly check that realistic amounts of primordial non- 
Gaussianity have no significant effect on the density profile 
of this saddle-point. This ensures that shell-crossing ap- 
pears for (almost) the same nonlinear density 5+ ^ 200 
(see Valageas 2009b), so that the large- mass tail of the halo 
mass function can be derived provided halos are defined 
by a nonlinear density threshold that is below this upper 
bound (which is indeed the case). Although this procedure 
only gives the large-mass tail, we have proposed a simple 
change of variable, applied to the mass function fitted to 
Gaussian numerical simulations, that obeys this large-mass 
asymptotic while keeping the mass function normalized to 
unity. If one uses the Press-Schechter mass function this 
gives back the result of Matarrese et al. (2000) but we argue 
that this procedure is somewhat more natural if one wishes 
to recover a more accurate mass function in the Gaussian 
case. We also checked that this agrees well with results from 
non-Gaussian numerical simulations. 

Next, we applied this method to the two-point corre- 
lation of massive halos, following the approach of Kaiser 
(1984). As in Valageas (2009b) we take into account the 
displacement of halo pairs under their mutual gravitational 
attraction. This gives the real-space halo correlation ^M(a;), 
whence the real-space bias 6m (a;). Since this approach does 
not assume that the halo correlation is small, the nonlin- 
ear formula it yields can be used to check whether the 
"linearized" form (where one only keeps the linear term 
over the matter correlation £,{x)) is a good approximation 
in the regime of interest. As expected, we find that the 
correction A6M(a;, /nl) to the Gaussian bias grows with 
6Af(a;,0) and with scale, roughly as /NL6M(a;, 0)x^, up to 
X ~ 100/i~^Mpc. Beyond this scale, the baryon acoustic os- 
cillation and the fact that the halo and matter correlations 
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do not change sign at the same point lead to strong oscilla- 
tions and divergent spikes for hM{x). This means that for 
X > 100/i~^Mpc the bias is no longer a useful quantity and 
one should directly work with the halo and matter corre- 
lations. In agreement with Desjacques (2008) we find that 
the two-point correlation of massive halos, that have a large 
bias, strongly amplifies the baryon acoustic oscillation. In 
addition we also obtain the modifications associated with 
primordial non-Gaussianity. The baryon oscillation remains 
strongly amplified, with a small shift, but somewhat less so 
for positive /nl- 

Finally, we used the "linearized" form of the halo two- 
point correlation to derive the halo power spectrum and the 
halo bias in Fourier space. We also give a simple recipe that 
ensures that the halo power spectrum always remains posi- 
tive (this only differs from the direct prediction by terms of 
order /^^ ^^id higher). We obtain a good agreement with 
numerical simulations, without introducing any free param- 
eter. Moreover, the two formulae described above allow one 
to estimate the range over which linear approximations over 
/nl are sufficient. Thus, we find that terms of order 
start playing a role at low k {k < 0.01/iMpc~^) for large 
negative /nl (/nl < —100), where the direct formula would 
give a negative power spectrum. 

These results, which do not involve free parameters 
(except for the mass function, if one requires its full 
shape, where one needs the fit to numerical simulations 
for Gaussian initial conditions) should be useful to con- 
strain primordial non-Gaussianities from observations of 
large-scale structures. In particular, as discussed above for 
Eg. ((95|) . the method presented in this article is quite gen- 
eral and can be applied to a large class of models. Moreover, 
since it provides results in both real space and Fourier space 
(i.e. the halo two-point correlation and power spectrum) it 
gives a complete and consistent description of halo clus- 
tering. As for previous approaches, the most reliable use of 
these models to constrain cosmology is to take advantage of 
the specific shape of the dependence on mass (for the mass 
function) or scale (for the bias) brought by primordial non- 
Gaussianity to constrain /nl (rather than the change in the 
amplitude at a given mass or scale). 
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